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PREDICTION  OF  PERFORMANCE  AND  DESIGN 

Spyros  A.  Kinnas 

Ocean  Engineering  Group,  Department  of  Civil  Engineering 
The  University  of  Texas  at  Austin,  Austin,  TX  78712,  USA 
http://cavity.ce.utexas.edu,  email:  kinnas@mail.utexas.edu 


ABSTRACT 

Recent  numerical  techniques  for  the  prediction  of  cavitating  flows,  in  linear  and  non-linear  theories,  are  applied 
on  super- cavitating  2-D,  3-D  hydrofoils  and  propellers.  Some  of  these  techniques,  when  incorporated  within  a 
non-linear  optimization  algorithm,  can  lead  to  efficient  super- cavitating  hydrofoil  or  propeller  designs.  This 
lecture  will  address  2-D  supercavitating  hydrofoils. 


1 INTRODUCTION 

High-speed  hydrofoil  or  propeller  applications  can  benefit  considerably,  in  terms  of  efficiency,  by  operating 
under  super-cavitating,  ventilating,  or  surface-piercing  conditions,  as  shown  in  Fig.  1.  A photograph  from  a 
super-cavitating  hydrofoil  experiment  that  was  performed  at  MIT’s  variable  pressure  Marine  Hydrodynamics 
Water  Tunnel  is  shown  in  Fig.  2.  Photographs  of  a modern  surface-piercing  propeller  when  out  and  in  the 
water  are  shown  in  Fig.  3. 

Cavitating  or  free-streamline  flows  were  first  addressed  in  nonlinear  theory  via  the  hodograph  technique 
as  introduced  by  Helmholtz,  Kirchoff  and  Levi-Civita  (Birkhoff  & Zarantonello  1957)1.  The  cavity  surface  in 
steady  flow  was  taken  as  a steamline  with  constant  pressure  (thus,  constant  velocity).  The  first  problems  to  be 
addressed  involved  flows  around  bluff  bodies  at  zero  cavitation  number. 

The  formulation  of  the  cavitating  flow  around  bodies  at  non-zero  cavitation  numbers  created  a lot  of 
diversity  on  the  cavity  termination  models.  Some  of  the  first  known  models  were  the  Riabouchinsky  end-plate 
model  and  the  re-entrant  jet  model.  Open  wake  models  were  introduced  in  an  attempt  to  model  the  viscous 
cavity  wake  (Fabula  1962,  Tulin  1964,  Yamaguchi  & Kato  1983,  Rowe  & Blottiaux  1993,  Kato  1994).  Some  of 
the  existing  cavity  termination  models  were  already  covered  in  previous  lectures  in  this  course. 

Due  to  the  difficulty  of  the  hodograph  technique  in  dealing  with  general  body  shapes,  very  few  cases  have 
been  treated  analytically.  The  hodograph  technique  was  extended  numerically  to  treat  arbitrary  geometries 
(Wu  & Wang  1964)  and  later  applied  to  the  analysis  of  super-cavitating  hydrofoils  in  the  presence  of  a free 
surface  (Furuya  1975a).  This  method  however,  still  could  not  treat  general  shape  three  dimensional  geometries. 

The  linearized  cavity  theory  was  introduced  by  (Tulin  1953)  and  became  quickly  very  popular,  as  proven 
by  the  vast  amount  of  publications2  which  made  use  of  it. 

Unfortunately,  the  linearized  theory  for  partially  cavitating  hydrofoils  predicts  that  by  increasing  the  thick- 
ness of  a hydrofoil,  the  extent  and  size  of  the  cavity,  for  constant  flow  conditions,  also  increases.  This  contradicts 
the  fact  that  thicker  hydrofoils  have  larger  leading  edge  radii  which  are  known  to  delay  cavitation  inception 
and  therefore  are  expected  to  develop  smaller  cavities.  In  addition,  it  is  well  known  that  linear  theory  grossly 
overpredicts  the  cavity  extent  and  volume,  especially  in  the  case  of  partial  cavitation. 

(Tulin  & Hsu  1980)  developed  the  short  cavity  theory  by  considering  the  cavitating  flow  as  a small  per- 
turbation on  the  nonlinear  fully  wetted  flow.  Thus,  the  nonlinear  foil  thickness  effects  were  included  in  this 

1The  list  of  references  is  located  at  the  end  of  the  second  lecture  of  Prof.  Kinnas. 

2An  extended  list  of  which  may  be  found  in  (Tulin  & Hsu  1980)  or  (Kinnas  1991). 


Paper  presented  at  the  RTOA  VT  Lecture  Series  on  "Supercavitating  Flows  ”,  held  at  the  von  Kdrmdn 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010. 
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Super-cavitating  Hydrofoil  Surface-piercing  Hydrofoil 


Figure  1:  Schematic  of  different  types  of  high-speed  hydrofoils. 


Figure  2:  Photo  of  a supercavitating  hydrofoil  experiment  inside  MIT’s  Marine  Hydrodynamics  Water  Tunnel, 
l/c&  3,  from  Kinnas  & Mazel,  1993. 
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formulation.  This  method  predicted  that  by  increasing  the  thickness  of  a partially  cavitating  hydrofoil,  the  size 
of  the  cavity  was  reduced  substantially  for  fixed  flow  conditions. 

A nonlinear  numerical  method  was  employed  to  analyze  cavitating  hydrofoils  by  using  surface  vorticity 
techniques  and  by  applying  the  exact  boundary  conditions  on  the  cavity  and  on  the  foil  (Uhlman  1987,  Uhlman 
1989).  An  end-plate  cavity  termination  model  was  implemented.  A reduction  in  the  size  of  the  cavity  as  the  foil 
thickness  increased  was  predicted,  but  not  as  drastic  as  that  predicted  in  (Tulin  & Hsu  1980).  A surface  vorticity 
technique  to  deal  with  thick  foil  sections  which  employed  an  open  cavity  model  was  developed  in  (Yamaguchi 
& Kato  1983).  Similar  boundary  element  method  techniques  were  developed  by  (Lemonnier  & Rowe  1988)  and 
by  (Rowe  & Blottiaux  1993).  Potential  based  boundary  element  methods  were  finally  applied  by  (Kinnas  & 
Fine  1991b,  Kinnas  & Fine  1993)  and  by  (Lee  et  al  1992). 

Three-dimensional  flow  effects  around  cavitating  finite  span  hydrofoils  were  treated  first  in  strip-theory  via 
matching  with  an  inner  two-dimensional  solution  within  either  linear  (Nishiyama  1970,  Leehey  1971,  Uhlman 
1978,  Van  Houten  1982)  or  non-linear  theory  (Furuya  1975b). 

The  complete  three-dimensional  super-cavitating  hydrofoil  problem  was  first  treated  in  linear  theory  via  a 
numerical  lifting  surface  approach  based  on  the  pressure  source  and  doublet  technique  (Widnall  1966)  and  later 
via  a vortex  and  source  lattice  technique  (Jiang  & Leehey  1977).  In  the  latter  work,  an  iterative  scheme  was 
introduced  which  determined  the  extent  of  the  cavity  by  requiring  the  pressure  distribution  on  the  cavity  to  be 
constant  along  the  span  (in  addition  to  being  constant  along  the  chord).  A variational  approach  for  determining 
the  cavity  planform  was  introduced  in  (Achkinadze  & Fridman  1994). 

Numerical  boundary  element  methods  within  non-linear  cavity  theory  were  naturally  extended  to  treat 
super-cavitating  3-D  hydrofoils  (Pellone  & Rowe  1981)  and  3-D  hydrofoils  with  partial  cavities  (Kinnas  & Fine 
1993)  or  cavities  with  mixed  (partial  and  super-cavities)  planforms  (Fine  & Kinnas  1993a).  Similar  methods 
were  also  developed  by  (Kim  et  al  1994,  Pellone  & Peallat  1995). 

The  first  effort  to  analyze  the  complete  three  dimensional  unsteady  flow  around  a cavitating  propeller 
subject  to  a spatially  non-uniform  inflow  was  presented  in  (Lee  1979,  Lee  1981,  Breslin  et  al  1982).  A source 
and  vortex  lattice  lifting  surface  scheme  was  employed  and  the  unsteady  three  dimensional  linearized  boundary 
conditions  were  applied  on  the  cavity.  The  cavity  planform  was  determined  at  each  blade  strip  and  each  time 
step  (i.e.,  blade  angle)  by  searching  for  the  cavity  length  which  would  produce  the  desired  vapor  pressure  inside 
the  cavity.  The  effect  of  the  other  strips  was  accounted  for  in  an  iterative  sense  by  “sweeping”  along  the  spanwise 
direction  of  the  blade  back  and  forth  until  the  cavity  shape  converged.  Similar  methods  were  presented  more 
recently  by  (Ishii  1992,  Szantyr  1994,  Kudo  & Ukon  1994). 

Unfortunately,  all  these  3-D  methods  are  hampered  by  the  inherent  inability  of  linear  cavity  theory  to  predict 
the  correct  effect  of  blade  thickness  on  cavity  shape,  as  already  mentioned.  This  deficiency  was  corrected  in  two 
dimensions  in  (Kinnas  1985,  Kinnas  1991),  where  the  leading  edge  correction  was  introduced  in  the  linearized 
dynamic  boundary  condition  on  the  cavity.  The  leading  edge  correction  was  subsequently  applied  to  the  three 
dimensional  propeller  solution  (Kerwin  et  al  1986,  Kinnas  1992b). 

Non-linear  methods  based  on  an  assumed  semi-elliptic  cavity  sectional  shape  have  also  been  applied  (Stern 
& Yorus  1983)  and  (Van  Gent  1994).  Non-linear  potential-based  boundary  element  methods  were  finally  applied 
to  cavitating  propellers  in  non-uniform  flows  by  (Kinnas  & Fine  1992,  Fine  & Kinnas  1993b),  and  more  recently 
by  (Kim  & Lee  1996). 

The  inviscid  cavity  flow  method  was  coupled  with  a boundary  layer  solver  in  the  case  of  partial  and  super- 
cavitating  2-D  hydrofoils  by  (Kinnas  et  al  1994).  This  allowed  for  the  inclusion  of  the  viscous  boundary  layer 
in  the  wake  of  the  cavity  and  for  determining  the  cavity  detachment  point  based  on  the  viscous  flow  upstream 
of  the  cavity. 

Reynolds- Averaged  Navier-Stokes  solvers  have  also  been  applied  in  the  case  of  the  prediction  of  attached 
sheet  cavitation  on  2-D  hydrodoils  (Kubota  et  al  1989,  Deshpande  et  al  1993).  An  overview  of  viscous  flow 
solvers  applied  to  cavitating  flows  may  be  found  in  (Kato  1996).  However,  these  methods  appear  to  be  best 
suited  for  the  prediction  of  cloud  and  or  detached  cavitation. 

In  these  two  lectures,  linear  and  nonlinear  methods  for  the  prediction  of  super-cavitation  on  hydrofoils 
and  propellers  will  be  summarized,  and  some  comparisons  with  experiments  will  be  presented.  Non-linear 
optimization  techniques,  applicable  to  the  design  of  2-D  super-cavitating  sections  and  super-cavitating  propeller 
blades,  will  be  also  presented. 
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Dynamic  Boundary  Condition 

Cavity  Detachment  (constant  pressure  condition)  Cavity  Closure  Condition 


Figure  4:  Formulation  of  the  inviscid  cavity  flow  problem 


2 2-D  HYDROFOIL 

2.1  Formulation 

Consider  the  2-D  super-cavitating  hydrofoil3,  as  shown  in  Fig.  4.  Assuming  inviscid  and  irrotational  flow, 
the  governing  equation  everywhere  inside  the  fluid  region  is  given  by4: 


v'U  = o (i) 

where  <j>  is  the  the  perturbation  potential  defined  from: 

q = Uoc  + V<t>  (2) 

where  q is  the  total  velocity  vector  in  the  flow.  In  order  to  uniquely  determine  <j>,  the  following  boundary 
conditions  are  imposed: 

• On  the  wetted  foil  surface,  the  following  kinematic  boundary  condition  is  applied,  which  requires  the  fluid 
flow  to  be  tangent  to  the  surface  of  the  foil.  Therefore, 


dcj> 

dn 


= -Uc 


n 


(3) 


where  n is  the  surface  unit  normal  vector. 

• At  infinity  the  perturbation  velocities  should  go  to  zero. 

V<{>  -»  0 (4) 

• The  dynamic  boundary  condition  specifies  constant  pressure  on  the  cavity,  or  (via  Bernoulli  equation) 
constant  cavity  velocity  qc: 


qc  = tWl  +cr  (5) 

where  the  cavitation  number,  a,  is  defined  as: 


a = 


Poo  Pc 

£772 

2uoo 


(6) 


Poo  is  the  pressure  corresponding  to  a point  in  the  free-stream  and  pc  is  the  pressure  inside  the  cavity. 

3The  application  of  the  methods  to  partially  cavitating  hydrofoils  is  straight-forward  and  the  reader  can  find  more  details  in 
(Kinnas  1998). 

4The  cavity  is  assumed  to  detach  at  a known  location  on  the  foil.  A criterion  for  determining  this  location  will  be  discussed  in 
a later  section. 
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• The  following  conditions  at  the  trailing  edge  of  the  cavity: 

1.  The  cavity  closes  at  its  trailing  edge.  The  inclusion  of  the  viscous  wake  downstream  of  the  cavity 
will  be  addressed  in  a later  section. 

2.  Near  the  trailing  edge  of  the  cavity,  the  pressure  recovery  termination  model  renders  for  the  cavity 
velocity,  qtr,  over  a transition  region  Ac: 

qtr  = UooVl  + v [1  - /(a 0]  (7) 

where  f(x)  is  an  algebraic  function  defined  in  (Kinnas  & Fine  1993). 

The  problem  of  finding  the  cavitation  number  for  given  cavity  extent,  l,  will  be  addressed  first  in  the  context 
of  linear  and  nonlinear  theories.  A method  for  determining  the  cavity  extent  for  given  cavitation  number  will 
be  described  in  a later  section. 

2.2  Linear  theories 

In  this  section,  the  linearized  cavitating  hydrofoil  problem  is  formulated  in  terms  of  singular  integral  equa- 
tions with  respect  to  unknown  vorticity  and  source  distributions.  These  integral  equations  are  inverted  analyt- 
ically in  the  case  of  general  shape  supercavitating  hydrofoils.  The  cavitation  number,  the  vorticity  and  source 
distributions  are  expressed  in  terms  of  integrals  of  quantities  which  depend  only  on  the  foil  geometry  and  the 
cavity  length.  The  procedure  is  summarized  here.  The  details  are  given  in  (Kinnas  1992a)  and  (Kinnas  & Fine 
1991a). 

We  define  as  u and  v the  perturbation  velocities  tangent  and  normal  to  the  direction  of  the  incoming  flow 
respectively,  as  shown  in  Figure  5.  In  the  context  of  the  linearized  cavity  theory  the  boundary  conditions  of  the 
corresponding  Hilbert  problem  are: 

The  kinematic  boundary  conditions: 
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U Uqq 


drjt  _ 
dx  ’ 


v+  = U00 


drju  . 
dx 


The  dynamic  boundary  conditions: 


„+  _ 


a 


0 < x < 1 , y = 0 
0 < x < Id  , y = 0+ 


Uoo  ‘f  Id  x l , y — 0 


a , 


u = — ?7°o  ; 1 < x < / , y = 0" 


(8) 

(9) 


(10) 

(11) 


where  r/i(x)  and  r}u(x)  are  the  ordinates  of  the  lower  and  upper  hydrofoil  surface,  respectively,  as  shown  in 
Figure  5.  These  boundary  conditions  can  be  expressed  in  terms  of  vorticity  and  source  distributions  5 y( x ) and 
q(x),  respectively,  located  on  the  slit  x € [0,/]. 


v-  = i , 1 f1 

2 27r  Jo  £-  x 

u+  = 

2 2ttJo 

u-  = 7 1 

2 2ttJo 

By  using  equations  (10),  (11),  (13),  and  (14)  it  can  easily  be  shown  that: 


j(x)  = 0;  1 < x < l 

Finally,  and  with  the  use  of  the  definitions: 


y(ar)  = 


7(g) 

crUoo 


and  q(x)  = 


gfo) 

wU oo 


the  complete  boundary  value  problem  becomes: 


1.  Kinematic  Boundary  Conditions 


(12) 

(13) 

(14) 

(15) 

(16) 


0 

1 


2 + ^f  = 0«W  5 0 < * < lD  , y = 0+ 


2.  Dynamic  Boundary  Condition 


1 

2 


1 fq(m 

2w  J £ — x 


o 


1 

2 


Id  <r  <1  , y = 0+ 


(17) 

(18) 


(19) 


3.  Kutta  Condition 0 


7(1)  = 0 (20) 

5With  -p  designating  the  Cauchy  principal  value  of  the  integral. 

®The  application  of  a Kutta  condition  may  seem  unnecessary  due  to  the  requirement  of  a finite  pressure,  thus  velocity,  at  the 
foil  trailing  edge.  Nevertheless,  this  condition  is  still  required  when  inverting  the  integral  equations. 
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4.  Cavity  Closure  Condition  ' 


where 


0? 


i 


1 drj, 
a dx 


/ q(x)dx  = 0 

(21) 

J 

0 

, r.*  1 dr)u 

■ and  = — 

“ a dx 

(22) 

2.2.1  Inversion  of  the  Integral  Equations  - Id  = 0 

In  the  case  of  leading  edge  detachment  (Id  = 0),  the  singular  integral  equations  of  Cauchy  type,  (17)  and 
(19),  can  be  inverted  to  produce  expressions  for  the  unknown  a,  y(x)  and  q(x)  in  terms  of  the  cavity  length  I 
and  the  lower  hydrofoil  surface  r/i(x),  as  follows  (Kinnas  1992a): 

First,  equation  (19)  is  inverted  with  respect  to  the  unknown  q(x)  (Muskhelishvili  1946)  to  produce: 


q(x)  = - 


+ 1 [JEL 1 AZZim 

7tV  l-XJo  V £ 


(23) 


where  use  of  equation  (15)  has  been  made.  Notice  that  the  expression  (23)  corresponds  to  the  unique  solution 
to  (19)  which  behaves  like  1 j\/l  — x at  the  trailing  edge  of  the  cavity  ( Wu’s  singularity  (Wu  1957)  ).  By 
substituting  equation  (23)  in  (17)  and  by  using  the  substitutions: 


z = 


we  arrive  at  the  following  singular  integral  equation  of  Cauchy  type  for  7: 


J_  fl  7 (V)dy  _ z ®*t(z) 

27 rJl  (1  + i)2 *)(z  - rf)  4(1  + z2)  2(1  + z2) 

Inversion  of  equation  (25)  with  respect  to  the  variable  7(77) /(I  + rj2)  renders  finally: 


(24) 


(25) 


7 (z)  = ~ 


(1  +z2)  ft^~z 


7T 


7' 


V 


f -erw 


-dij 


(26) 


t — rf  (l  + r)2)(z  - rj) 

Notice  that  7 (z)  in  equation  (26)  is  the  unique  solution  to  (25)  which  satisfies  the  Kutta  condition  (20)  at  z = t. 

The  cavitation  number  o is  determined  by  satisfying  the  cavity  closure  condition  (21).  First,  by  substituting 
equation  (23)  in  (21)  and  by  using  equation  (15)  we  can  get  (Kinnas  1992a): 


7 m = 0, 


(27) 


and  by  substituting  equation  (26)  we  can  get  the  following  general  expression  for  a (Kinnas  1992a): 


4y^2r4  f*  I rj  y/r2  + 1 + rjy/r2  — 1 

7 r(r2  + 1)  J0  y t - rj  (1  + rj2)2 


dm 

dx 


drj 


(28) 


where:  r = y/1  + t2. 

The  source  distribution  can  be  derived  by  substituting  equation  (26)  into  equation  (23)  (Kinnas  1992a): 


q{z) 


-0i  (*)  + 


t + z (y/r2  — 1 — zy/r2  + 1) 


2y/2 r2 


1 + z2  It  + z I ~to  @^(c u)cku 

7T  V Z J0  V t ~ CU  (1  + IU2)(Z  + u) 


(29) 


rWe  apply  the  linearized  cavity  closure  condition  in  which  the  cavity  is  required  to  have  zero  thickness  at  its  trailing  edge.  The 

present  method  can  be  extended  to  treat  open  cavities  at  the  trailing  edge  with  the  openness  of  the  cavity,  possibly  supplied  from 
further  knowledge  of  the  viscous  wake  behind  the  cavity.  The  effects  of  viscosity  will  be  addressed  at  a later  section. 
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for  z < t,  and: 


q(z) 


t + z (yV2  — 1 — zsjr2  + 1) 


z 

1 + z2  lt  + z 


2\/2r 


7T 


j*  CJ 


uj  &i(u)du 


c o (1  + CU2)(z  + uj) 


— t (, z\/r 2 + 1 + i Jr'2  — 1) 


2v/2 r 


1 + z 2 /z-t 


7T 


/* 

JO 


uj  &i(u)doj 


t — UJ  (1  + uj2)(z  — uj) 


(30) 


for  z > t. 


2.2.2  Inversion  of  the  Integral  Equations  - Id  > 0 

In  this  case,  as  shown  in  (Kinnas  & Fine  1991a),  equations  (26),  (29),  and  (30)  still  apply,  after  the  following 
substitutions  are  made: 


and 


replace  7 with  7 — 2 (u+w  — -)  ; 0 < x < Id 


(31) 


replace  &*  with  0;*  + F 


(32) 


with 


U + 1 

F(x)  = --  f U\ 7 ~ 2 <%  (33) 

ttJ  — x 
0 

where  is  the  horizontal  perturbation  velocity  on  the  wetted  part  on  the  suction  side  of  the  foil,  divided  by 
aUoo-  The  value  of  u+w  is  determined  by  applying  the  kinematic  boundary  condition  on  the  upper  wetted  part 
of  the  hydrofoil,  equation  (9).  This  is  equivalent  to  requiring  that  the  value  of  q(x)  for  0 < x < Id  is  equal 
to  the  value  of  the  thickness  source  in  the  case  of  wetted  flow.  A rather  lengthy  formula  for  u+,,  as  well  an 
expression  for  the  modified  value  of  a are  given  in  (Kinnas  & Fine  1991a),  and  are  not  included  in  this  lecture. 


2.2.3  The  cavity  shape 

The  cavity  thickness  h(x),  which  also  includes  the  foil  thickness  as  shown  in  Figure  5,  is  determined,  within 
the  framework  of  linearized  theory,  by  integrating  the  equation: 

Uoo^  = ^ 

The  camber  of  the  cavity  in  the  wake,  c(x),  is  determined  by  integrating  the  following  equation: 

dc 

Uoo~r  = vw(+)  for  1 < x < l (35) 

ax 

where  vw(x)  is  the  normal  perturbation  velocity  in  the  wake,  given  as  follows: 

Vw(x)  = --^-[  ; l<x<l,  y = 0 (36) 

"T  Jo  X ~ S 

By  substituting  equation  (26)  in  (36)  we  can  get  vw  in  terms  of  the  hydrofoil  geometry  (Kinnas  & Fine 
1991a): 


21-10 


vw (z)  __  jt  + z (Vr2  - 1 - z\/r2  + 1) 

crlJoo  V z 4\/2 r2 

1+Z2  It  + Z f*  I U 0;*  (oj)dhj 

2?r  V z J0  ]/t  — u (1  +u2)(z+u) 

jz  — t {z\Jr2  + 1 + V r2  — 1) 

V ~ 4v/2r2 

1 + z2  lz  — t r*  I w 0* (w)du 

2?r  V ^ /0  yi— tu  (1  + w2)(;j  — (jS) 


(37) 


The  pressure  distribution  on  the  upper  and  lower  cavity  or  foil  surface  is  given,  in  the  context  of  linear 
theory,  as  follows: 


C+  = -2a 


1 _ J_  £ qQdt 

2 2tt  J £ — x 


C~  = -2cr 


7 _ J_  f q(0<% 

'2  2tt  J Z-x 


where  Cp  is  the  pressure  coefficient  defined  as: 


0 < x <1  , y = 0+ 


0 < x < l , y = 0 


^ P-Poo 

- urn 

2 uoo 


(38) 


(39) 


(40) 


2.2.4  Numerical  Integrations 

The  integrals  in  equations  (28),  (26),  (29)  and  (30)  are  computed  numerically  with  special  care  taken  at 
the  singularities  of  the  integrands.  We  first  define  the  transformation: 


rj  = t 


0 < p < t and  0 <9  < tt 


(41) 


Next,  we  expess  the  involved  integrals  in  terms  of  9 , thus  avoiding  the  square  root  singularities  of  the 
integrands  at  rj  = 0 and  rj  = t . The  numerical  integrations  are  then  performed  by  applying  Simpson’s  rule  with 
K uniform  intervals  in  9. 

To  compute  the  principal  value  of  the  singular  integral  in  equation  (26),  we  first  factor  out  the  involved 
singularity  as  follows  (Kinnas  1992a): 


/(*?)  ~ /(*) 
Z — Tf 


dr)  - Trf(z) 


(42) 


where  f{rj)  = ^7+^^  • Notice  that  the  integrand  in  the  integral  of  equation  (42)  is  not  singular  anymore,  and 
thus  the  integral  is  computed  numerically  by  applying  the  same  methodology  described  in  the  beginning  of  this 
section.  As  shown  in  (Kinnas  1992a),  Simpson’s  rule  produces  very  accurate  values  for  the  integrals  even  with 
five  uniform  intervals  ( K = 5). 


2.2.5  The  vortex/source-lattice  method 

A direct  numerical  method  must  be  applied  in  the  case  of  3-D  hydrofoils  or  propeller  blades,  as  in  (Lee 
1979,  Kinnas  et  al  1998a,  Kosal  1999).  In  this  case  the  involved  integrals  are  disretized  first,  then  the  boundary 
conditions  are  applied  at  some  appropriately  selected  control  points  (C.P.s),  and  finally,  the  resulting  system 
of  linear  equations  is  inverted.  The  integrals  in  equations  (17)  and  (19)  are  disretized  over  each  segment  i 
by  replacing  7 (£)d£  with  ji(Xpv.+1  — XpVi),  and  q{£)dt;  with  qi(XpS[+1  — XpSi).  Figure  6 shows  two  types  of 
arrangements  for  the  discrete  vortices  and  sources  and  the  correpsonding  control  points. 
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U 

inflow 


Figure  6:  Half-cosine  and  full-cosine  discretization  scheme  in  the  vortex/source-lattice  method.  The  latter 
allows  for  a better  representation  of  typical  camber  lines  for  super-cavitating  hydrofoils. 


Figure  7:  Super-cavitaitng  non-linear  theory.  Definition  of  cavity,  (wetted)  body,  and  wake  surfaces.  The  cavity 
sources  in  the  wake  are  also  shown  in  the  case  of  the  hybrid  scheme. 
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Figure  8:  Super-cavitating  hydrofoil.  Definition  of  main  parameters.  Panel  arrangement  on  the  cavity  and  foil 
shown  for  N = 80. 


2.3  Nonlinear  theories 

A potential  based  boundary  element  method  which  has  applied  has  been  applied  for  the  analysis  of  cavitating 
hydrofoils  in  nonlinear  cavity  theory  (Kinnas  & Fine  1991b,  Kinnas  & Fine  1993,  Fine  & Kinnas  1993a),  is 
summarized  here,  with  emphasis  on  super-cavitating  hydrofoils.  The  perturbation  potential  on  the  combined 
foil  and  cavity  surface  satisfies  the  following  integral  equation  (Green’s  third  identity): 


7T(f)p 


Y~\n  R-<j> 
on 


9 In  R 
dn 


dS 


din  R' 
dn 


dS 


(43) 


n is  a unit  vector  normal  to  the  foil  or  cavity  surface,  Swb  is  the  wetted  body  (foil)  surface,  Sc,  the  cavity 
surface,  and  Sw  is  the  trailing  wake  surface,  as  shown  in  Fig.  7.  R is  the  distance  between  a point,  P,  and  the 
point  of  integration  over  the  wetted  foil,  cavity,  or  wake  surface,  as  shown  in  Fig.  7. 

The  foil  and  cavity  surface  are  discretized  into  flat  panels,  as  shown  in  Fig.  8.  The  source  and  dipole 
strengths  are  assumed  constant  over  each  panel.  On  the  wetted  foil  surface,  the  source  strengths,  which  are 
proportional  to  d(f>/dn,  are  given  by  the  kinematic  boundary  condition,  equation  (3).  On  the  cavity,  the  dipole 
strengths,  which  are  proportional  to  <f>,  are  determined  from  the  application  of  the  dynamic  boundary  condition 
(5)  and  (7).  For  simplicity  A = 0 for  the  rest  of  this  paper;  the  complete  derivation  is  given  in  (Kinnas  & 
Fine  1993).  It  should  be  noted  that  the  length  of  the  transition  region  A has  been  found  to  affect  the  results 
only  locally.  In  the  case  of  A = 0 and  when  a large  number  of  panels  is  used  the  formation  of  a re-entrant  is 
observed(Krishnaswamy  1999). 


qc  = Uoo  • s + = t/ooVTT?  (44) 

where  s is  the  arclength  along  the  cavity  surface  (measured  from  the  cavity  detachment  point),  and  s is  the  unit 
vector  tangent  to  the  cavity  surface.  Integration  of  equation  (44)  renders  the  potential  on  the  cavity  surface: 


<f>(s)  = 4>{ 0)  - Uoo  • s + sUoo-s/l  + n (45) 

Where  <j>{ 0)  is  the  potential  at  the  leading  edge  of  the  cavity.  In  the  numerical  scheme  4>{ 0)  is  expressed  in  terms 
of  the  (unknown)  potentials  on  the  wetted  part  of  the  foil  in  front  of  the  cavity. 

The  potential  jump,  A <j>m,  in  the  wake  is  determined  via  the  following  condition: 
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A(f>w  — 4>CTE  &CTE  (46) 

where  4>qTE  and  4>cte  are  the  potentials  at  the  upper  and  lower  cavity  trailing  edge  panels,  respectively.  In 
two  dimensions  the  effect  of  the  trailing  wake  surface  is  equivalent  to  the  effect  of  a concentrated  vortex  at  the 
cavity  trailing  edge  with  strength  equal  to  A<f>w. 

The  integral  equation  (43)  is  applied  at  the  panel  mid-points  together  with  equations  (46).  The  resulting 
linear  system  of  equations  is  inverted  in  order  to  provide  the  unknowns:  (a)  <j>  on  the  wetted  foil,  (b)  d<j>/dn 
on  the  cavity,  and  (c)  the  corresponding  cavitation  number  a.  The  cavity  shape  is  determined  in  an  iterative 
manner.  In  the  first  iteration  the  panels  representing  the  cavity  are  placed  on  the  foil  surface  directly  under  the 
cavity.  In  subsequent  iterations  the  cavity  shape  is  updated  by  an  amount  h(s)  (applied  normal  to  the  cavity 
surface)  which  is  determined  from  integrating  the  following  ordinary  differential  equation8 : 


Qc 


dh 

ds 


— Uoo  n T 


d(f> 

dn 


In  addition,  the  cavity  closure  condition  is  enforced  by  the  following  equation: 


(47) 


'*(-)  = 7 / 


■n  + 


dcf> 

dn 


ds  = 0 


(48) 


sl  is  the  total  arclength  along  the  cavity  surface. 

The  predicted  cavity  shapes  and  cavitation  numbers  have  been  found  to  converge  quickly  with  number  of 
iterations,  especially  in  the  case  of  super-cavitating  hydrofoils. 

This  particular  feature  of  the  presented  method  makes  it  very  attractive  for  3-D  and/or  unsteady  flow 
applications,  where  carrying  more  than  one  iterations  would  increase  the  computation  time  substantially.  In 
fact  the  first  iteration  in  the  iterative  scheme  is  determined  by  using  the  following  hybrid  scheme. 


2.3.1  The  hybrid  scheme 

In  the  case  of  3-D  hydrofoils  and  propeller  blades,  the  following  hybrid  (combination  of  panel  and  source 
lattice  method)  scheme  has  been  developed  (Fine  & Kinnas  1993a).  In  this  scheme  the  panels  representing  the 
cavity  are  placed  either  on  the  foil  surface  underneath  the  cavity,  Scb,  or  on  the  approximate  wake  surface, 
Sew,  as  shown  in  Fig.  7.  Equation  (43)  then  becomes(Fine  & Kinnas  1993a): 


7T  < Pp  — 


+ / qw  In  RdS 

J Sew 


d(j>,  „ , d In  R 


SssWn^11  * dn 

f A cbw^^dS 

JScwUSw  °n 

for  P € SE 


dS 


(49) 


and 


27 T(f)p  = 7tA  4>w  + 

+ / qw  In  RdS  — 

J Sew 


LM 
L 


■In  R — (f> 


dlnR 


dn 


dS 


Scw^Sw 
for  P 6 Sew 


dlnR 

A<j)w  dS 


(50) 


where  Sb  = Swb  U Scb  is  the  surface  of  the  whole  foil,  and  qw  is  the  cavity  source  in  the  wake  surface  given 
as: 


d<t>+  d(j>  dhw 
dn  dn  ^ c ds 


(51) 


8As  shown  in  (Kinnas  & Fine  1993)  this  equation  is  equivalent  to  the  kinematic  boundary  condition  on  the  cavity  surface. 
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Figure  9:  Predicted  cavity  trailing  edge  thickness  and  cavity  shape  for  various  cavity  lengths  at  cr  = 1.097 
(corresponding  to  / = 0.4).  From  Kinnas  and  Fine,  1993. 


where  hw  is  the  cavity  height  measured  normal  to  Sew,  as  shown  in  Fig.  7.  The  cavity  height  h is  measured 
normal  to  the  foil  surface  for  the  part  of  the  cavity  that  overlaps  the  foil,  as  also  shown  in  Fig.  7. 

The  major  advantage  of  the  hybrid  scheme  is  that  it  can  apply  to  all,  wetted,  partially,  and  super- cavitating 
flows  alike  by  utilizing  the  same  panel  discretization.  In  fact  this  scheme,  as  described  in  later  sections,  is 
applied  on  super-cavitating  3-D  hydrofoils  and  propellers.  In  addition,  this  scheme  has  been  found  to  predict 
the  expected  non-linear  effect  of  foil  thickness  on  the  cavity  shape  in  the  case  of  partial  cavitation  (Kinnas  & 
Fine  1993). 

2.4  Cavity  extent  for  given  cavitation  number 

In  the  previous  sections  the  cavity  length  was  assumed  to  be  known.  In  the  case  of  linear  theory,  equation 
28  can  be  inverted  with  respect  to  l for  given  a.  In  the  case  of  nonlinear  theories  the  method  can  still  be  applied 
for  the  given  <j  and  various  “trial”  cavity  lengths.  In  this  case  though,  due  to  the  fact  that  the  cavitation  number 
is  given  (instead  of  being  determined),  the  cavity  closure  condition,  equation  (48),  or  its  equivalent  hw(l)  = 0 
in  the  case  of  the  hybrid  scheme,  will  not  be  satisfied. 

5(1-, a)  = hw(l ) 7^  0 (52) 

This  is  shown  for  a partially  cavitating  hydrofoil  in  Fig.  9 where  the  predicted  5 and  the  cavity  shapes  for 
fixed  a and  different  values  of  cavity  length  are  shown.  An  iterative  scheme  has  been  developed  (Kinnas  & Fine 
1993)  for  determining  the  cavity  extent  for  given  a by  searching  for  the  cavity  length,  l,  for  which  the  cavity 
closes  (within  a specified  tolerance)  at  its  trailing  edge9: 

6(1- a)  = 0 (53) 

A comparison  of  predictions  from  applying  the  linear,  the  fully  non-linear,  and  the  hybrid  cavity  models  to 
a super-cavitating  hydrofoil  is  shown  in  Figure  10.  All  theories,  including  the  conventional  linear  theory,  appear 


9 An  open  cavity  model  can  be  readily  implemented  within  this  method  by  requiring  the  specified  thickness  at  the  cavity  end. 
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Figure  10:  Comparison  of  predicted  cavity  shapes  from  linear  (dashed),  non-linear  (solid)  and  hybrid  (solid 
with  dots)  theories  for  NACA16004  at  a = 6°.  (dotted)  linear;  a/a  = 0.44.  From  Fine  and  Kinnas  (1993a). 


to  predict  the  cavity  shape  within  satisfactory  accuracy.  Note  that  linear  theory  overpredicts  the  slope  of  the 
cavity  at  the  leading  edge. 

Finally,  a comparison  of  predicted  versus  measured  velocity  profiles  in  the  vicinity  of  a 2-D  supercavity 
is  shown  in  Figure  11.  In  the  prediction  the  wall  effects  have  been  included  by  using  a sufficient  number  of 
multiple  images  with  respect  to  the  horizontal  tunnel  walls  (Kinnas  & Mazel  1992). 


2.5  Effects  of  viscosity 

Consider  a super-cavitating  hydrofoil  in  viscous  flow,  as  shown  in  Figure  12.  We  first  solve  the  inviscid 
cavity  flow  in  non-linear  theory  (Kinnas  & Fine  1991b).  We  then  apply  the  boundary  layer  equations  on  the 
compound  foil,  confined  by  the  pressure  side  of  the  hydrofoil  and  the  cavity  boundary.  A zero  friction  condition 
is  applied  everywhere  on  the  cavity,  as  shown  in  Figure  12.  The  effects  of  viscosity  may  be  included  by  including 
the  following  term  on  the  right-hand  side  of  equation  43: 


+ f a In  Rds  (54) 

J SwbUSc^Sw 

where  the  blowing  sources,  a,  are  defined  from  (Hufford  et  al  1994): 


a = 


dUe6 * 
ds 


(55) 


where  Ue  is  the  edge  velocity  (in  this  case  the  same  as  the  velocity  predicted  by  the  panel  method  when  applied 
on  the  modified  body  that  includes  the  displacement  thickness)  and  5*  is  the  displacement  thickess,  which  is 
determined  from  the  boundary  layer  equations  (Drela  1989). 

Results  from  applying  this  method  are  shown  in  Figure  13.  The  resulting  boundary  layer  displacement 
thickness  is  shown  at  the  top  of  the  figure.  The  pressure  distributions  in  inviscid  and  viscous  flow  are  also 
shown.  Notice  that  viscosity  has  a very  small  effect  on  the  pressure  distribution.  In  other  words,  for  super- 
cavitation, the  cavitation  number  in  viscous  flow  for  given  cavity  length,  l,  is  practically  identical  to  that  in 
inviscid  flow.  It  should  be  noted  that  in  the  case  of  partial  cavitation  the  effects  of  viscosity  are  much  stronger 
than  in  the  case  of  supercavitation  (Kinnas  et  al  1994).  The  friction  coefficient,  Cf,  on  the  pressure  side  of  the 
foil  and  cavity,  is  also  shown  at  the  bottom  part  of  Figure  13. 

The  lift  and  drag  on  the  hydrofoil  are  evaluated  by  integrating  the  pressure  forces  acting  on  all  sides  of  the 
hydrofoil  (the  constant  cavity  pressure  is  applied  on  the  cavitating  sides  of  the  hydrofoil)  as  well  as  the  frictional 
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Figure  11:  Numerical  vs.  experimental  velocities  as  a function  of  the  vertical  distance  from  the  cavity  boundary 
at  different  chordwise  locations.  From  Kinnas  and  Mazel  (1993). 


Figure  12:  Super-cavitating  hydrofoil  with  its  boundary  displacement  thickness. 
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N 

o V/c2 

cL 

cD 

100 

0.145  0.365 

0.282 

0.0219 

160 

0.146  0.364 

0.287 

0.0223 

200 

0.146  0.363 

0.292 

0.0231 

Table  1:  Convergence  of  viscous  cavity  solution  (a,  V/c 2,  Cl  and  Co)  with  number  of  panels.  Super-cavitating 
hydrofoil;  r/c  = 0.045, /0/c  = 0.015,  p = 0.85,  a = 3°. 


Figure  15:  The  main  parameters  in  the  design  of  a super-cavitating  hydrofoil. 


forces  acting  on  the  wetted  side  of  the  hydrofoil.  The  convergence  of  the  cavity  solution  and  the  predicted  forces 
with  number  of  panels  is  given  on  Table  1. 

Finally  the  predicted  a,  Cl,  and  Co  vs.  I curves  are  shown  in  Figure  14,  for  a super-cavitating  section  in 
inviscid  flow  and  for  two  Reynolds  numbers3 * * * * * * 10.  The  super-cavitating  section  is  a combination  of  a NACA  4digit 
camber  form  (with  the  maximum  camber  fa  at  x = p)  and  a linear  thickness  form.  The  effect  of  viscosity  on 
lift  coefficient  is  shown  to  be  very  small. 


3 DESIGN  OF  SURER- CAVITATIN G SECTIONS 
3.1  Statement  of  the  problem 

We  must  determine  the  cavitating  hydrofoil  geometry  and  its  operating  condition  (angle  of  attack,  a),  which 
produces  the  minimum  drag,  D,  for  specified  design  requirements. 

The  parameters  that  define  the  geometry  of  a super-cavitating  hydrofoil,  also  shown  in  Figure  15,  are: 

• the  chord,  c 

• the  maximum  camber,  /0,  on  the  pressure  side 

• the  location  of  the  maximum  camber,  fp 

The  design  requirements  taken  into  consideration  in  this  paper  are  : 


10Based  on  l. 
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• Sectional  lift,  L0(N/m) 

• Uniform  forward  velocity  of  the  foil,  Uoo(m/s) 

• Cavitation  number,  <r0,  defined  as  usual: 


Poo  Pv 

w 


where  p is  the  fluid  density,  p^  is  the  ambient  pressure,  and  pv  is  the  vapor  pressure. 


• Minimum  section  modulus  of  the  foil,  Zmin 


• Acceptable  cavity  length,  l 

• Acceptable  cavity  volume,  V,  or  cavity  height,  h 


(56) 


The  condition  on  the  cavity  length  is  necessary  in  order  to  avoid  unstable  cavities,  usually  being  the  long 
partial  or  the  short  super-cavities.  The  cavity  volume/height  constraints  ensures  acceptable  positive  cavity 
thickness  (volume)  in  order  to  avoid  very  thin  cavities  (also  negative  thickness  cavities)  which  are  either  non- 
physical or  may  turn  into  harmful  bubble  cavitation. 


3.2  The  hydrodynamic  coefficients 

The  hydrofoil  lift,  L,  and  drag,  D,  acting  on  the  hydrofoil,  as  shown  in  Figure  15,  can  be  expressed  in  terms 
of  lift  and  drag  coefficients,  Cl  and  Cd,  respectively: 


L = \pUlcCL 

(57) 

D = \pUlcCD 

(58) 

The  drag  coefficient,  Cd,  may  be  decomposed  into  two  components. 

CD  = Cb  + CVD 

(59) 

where  ClD  is  the  inviscid  cavity  drag  coefficient,  and  C’[}  is  the  viscous  drag  coefficient. 

For  known  hydrofoil  geometry,  angle  of  attack,  a,  and  cavity  length,  l,  any  hydrodynamic  or  cavity  quantity, 
Q (Cl,  ClD,  cr,  Vjc 2, ....),  may  be  expressed  as  follows: 

Q = Q(a , nondimensional  foil  geometry,  1/c)  (60) 

Inviscid  analysis  methods  are  used  to  determine  Q.  The  analytical  formulas  of  (Hanaoka  1964)  could  have 
been  used.  These  expressions  though  are  based  on  linearized  cavity  theory,  which  has  been  found  to  overpredict 
the  cavity  length  and  size  substantially,  especially  in  the  case  of  partial  cavitation.  Instead,  in  this  work  the 
already  described  numerical  non-linear  cavity  analysis  method  is  utilized  (Kinnas  & Fine  1991b,  Kinnas  & Fine 
1993).  The  shape  of  the  cavity  surface  is  determined  in  non-linear  theory  iteratively,  with  the  use  of  a low-order 
potential  based  panel  method.  The  (inviscid)  forces  are  determined  by  integrating  the  pressures  along  the  foil 
surface. 

The  viscous  drag  is  determined  by  assuming  a uniform  friction  coefficient,  Cf,  over  the  wetted  part  of  the 
foil.  Cf  is  expressed  in  terms  of  the  Reynolds  number  (Re  = U^c/ v,  v : kinematic  viscosity)  via  the  ITTC 
formula  (Comstock  1967): 

^ _ 0-675 

Cf  (l  0910  Re  — 2) 2 ( } 

During  the  optimization  process,  the  chord  length  (also  the  Reynolds  number)  is  not  known,  thus  the  value 
of  Cf  is  updated  at  each  optimization  iteration. 

Since  only  the  lower  surface  of  the  foil  determines  the  hydrodynamics  of  supercavitating  flows,  the  thickness 
is  not  included  as  a parameter.  The  upper  surface  can  be  placed  anywhere  arbitrarily  inside  the  cavity.  Thus, 
when  computing  the  section  modulus  of  a foil,  the  upper  cavity  surface  is  considered  as  the  upper  surface  of  the 
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“compound”  foil.  It  is  reasonable  to  assume  that  the  cavity  always  starts  at  the  leading  edge  of  the  foil,  since 
the  supercavitating  sections  have  a sharp  leading  edge.  Furthermore,  we  deal  with  situations  where  the  cavity 
detaches  at  the  trailing  edge  of  the  foil  on  the  pressure  side.  For  the  case  of  supercavitating  sections,  in  addition 
to  specifying  the  minimum  section  modulus  of  the  compound  section,  the  minimum  allowable  cavity  height  at 
the  10%  of  the  chord  length  from  the  leading  edge  is  specified  via  an  inequality  constraint.  This  ensures  positive 
cavity  thickness,  as  well  as  sufficient  local  strength  at  the  sharp  leading  edge  of  the  foil. 

The  relevant  coefficients  are  expressed  as  follows: 

CL  = CL  (a,  /o  / c,  fp/c,  l /c) 

Cb  = Cb(a,fo/c,fP/c,l/c) 
a = a(a,fo/c,fp/c,l/c) 

V/c2  = V/c2(a,f0/c,fp/c,l/c) 
zf(?  = z/c3(a,  fo/c,  fp/c,l/c) 
hw/c  = hio/c(a,fo/c,fp/c,l/c) 

where  hio  is  the  cavity  height  at  the  10%  of  the  chord  length  from  the  leading  edge  of  the  foil.  Note  that  the 
section  modulus  is  a function  of  all  the  geometry  variables  including  the  cavity  length,  whereas  it  is  independent 
of  the  cavity  length  in  the  case  of  partial  cavitation. 

The  coefficients  are  computed,  as  in  the  case  of  partial  cavitation,  for  a number  of  foil  geometries  by  using 
the  analysis  method  of  (Fine  & Kinnas  1993a).  The  coefficients  are  again  expressed  by  using  second  degree 
polynomials  in  terms  of  all  the  geometric  parameters.  For  example, 


CL  = CLla  + CL2a2  + CLs^  + CLi 
+CL,  + CLe  +CLrl-  + CL,  (j) 

+CLta ^ + Cl10oc—  + CLlla—  + CLJ^ 

c c c c c 

+Cli3^+Cl^¥c+Cli5 

where, 

r i 

C “ ^1  - (c/02  + 1 - (c/02 

The  variable  transformation  from  l/c  to  V/c  is  done  knowing  that  Cl  and  Co  are  linear  in  V [c  for  the 
super-cavitating  flat  plate  in  linear  theory,  as  described  by  (Geurst  1960).  This  transformation  has  been  found 
to  improve  the  accuracy  of  the  interpolation  substantially.  The  effect  of  the  number  of  The  accuracy  of  the 
interpolation  for  Cd  may  be  seen  in  Figures  16. 

3.3  The  optimization  problem 

Consider  the  following  general  nonlinear  optimization  problem: 

minimize  /(x) 

subject  to  ^(x)  < 0 i = 1,2,..., m 

hi(x)  = 0 i = l,2,...,l  (64) 

where  /(x)  is  the  objective  function  defined  on  7 Zn.  x is  the  solution  vector  of  n components.  #i(x)  < 
0,  ■ ■ • , flw(x)  < 0 are  inequality  constraints  defined  on  fcn  and  hi (x)  = 0, . . . , hi(x)  = 0 are  equality  constraints 
also  defined  on  'Rn. 

The  solution  vector  is  defined  as  (n  = 5): 


(62) 

(63) 


x = [a,c,  fo,  fp,l]T 


(65) 
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The  constraints  are  (m  = 3,  l = 2) 


/ii(x) 

M*) 

51  (x) 

52  (x) 

53  (x) 


T(x)  -L0  = 0 
a(x)  - <70  = 0 


<0 


Zmin 

c3  c3 


< o 


/Ho(x)  | / /no A <() 


(66) 

(67) 

(68) 

(69) 

(70) 


3.4  The  numerical  optimization  method 

The  method  of  multipliers  is  applied  in  which,  first  each  of  the  inequality  contraints,  <?i(x)  < 0,  is  converted 
into  an  equality  constraint  with  the  introduction  of  the  new  variable,  s*: 

5;(x)  + 4=°  (71) 

Then  the  augmented  Langrangian  penelty  function  is  formed  (Mishima  & Kinnas  1996): 


F(x  i, ...,  «i, ...,  Ai, ...,  Ai,  ...)  = 

D(x)  + SA ihi(x.)  -)-  SA i[gi(x)  + sf] 

+Scj/i.f  (x)  + Sci[^(x)  + s?]2  (72) 

where  A*  and  A*  are  the  Langrange  multipliers,  and  c,  and  c,;  are  penalty  parameters.  The  parameters  of  the 
problem,  the  additional  parameters  («j),  the  Langrange  multipliers,  and  the  penalty  parameters  are  determined 
by  minimizing  function  F.  A technique  to  minimize  F is  described  in  great  detail  in  (Mishima  & Kinnas  1996). 

3.5  Effect  of  the  initial  guess  on  the  optimum  solution 

In  general,  a nonlinear  optimization  problem  is  initial  guess  dependent.  In  other  words,  if  there  exist  more 
than  one  local  minimum  to  the  objective  function,  any  one  of  these  minima  may  be  obtained  depending  on  the 
initial  condition.  Due  to  the  special  structure  of  the  problem,  no  multiple  solutions  were  found  for  the  range  of 
cavitation  numbers  that  we  studied.  This  uniqueness  is  attributed  partly  to  the  fact  that  the  involved  functions 
are  quadratic,  i.e.  not  wavy,  and  partly  to  the  fact  that  the  range  of  the  solutions  is  known,  so  that  a reasonable 
initial  guess  can  be  selected.  Figure  17  shows  that  several  different  initial  guesses  lead  to  the  same  solution. 


3.6  Results 

A NACA  4 digit  camber  form  and  the  Johnson  five-term  camber  form  are  used.  The  NACA  4digit  camber 
form  (Abbott  & Von  Doenhoff  1959)  has  two  parameters,  which  are  the  maximum  camber  /0/c  and  the  location 
of  the  maximum  camber  fp/c,  as  shown  in  Figure  15.  The  values  for  the  constraints  are: 


1.15 


hm 

c 


0.01 


( Zmin  \ 

v-^ ) 


7 x 10 


-5 


(73) 


The  method  is  applied  for  fixed  lift  (L0  = 30,000 N/m)  and  for  a range  of  cavitation  numbers  between 
0.15  and  0.7  (corresponding  to  approximate  speeds  of  70  and  30  knots,  respectively11)  for  partial  and  super- 
cavitating  conditions.  The  resulting  L[D  for  both  cases  are  shown  in  Figure  18.  As  expected,  for  “low”  a a 
supercavitating  section  has  larger  L[D  than  a partially  cavitating  section,  and  the  reverse  holds  for  larger  a. 
Optimum  partially  and  supercavitating  sections  are  also  shown  in  Figure  18. 

Finally,  we  show  in  Figure  19,  contour  plots  of  L/D  for  optimum  sections,  designed  by  the  present  method, 
over  a range  of  combinations  of  required  lift  and  cavitation  number.  This  graph  is  intended  to  help  the  designer 
decide  as  to  what  is  the  best  solution  (partially  or  supercavitating  section),  depending  on  the  requirements. 


nThe  cavitation  number  is  inversely  proportional  to  the  square  of  the  speed  of  the  hydrofoil,  for  fixed  ambient  pressure  and 
temperature  conditions. 


21-23 


3.7  Comparison  of  the  present  to  existing  design  methods 

The  designed  sections,  under  given  requirements,  are  compared  to  those  designed  from  other  methods. 
These  methods  are  based  on  either  the  Tulin  two-term  or  the  Johnson  five-term  sections,  given  next: 


V 

c 


(Tulin  two  — term) 


(74) 


/X\5/2  fX\Z  fX\7/2. 

-30912  (-J  + 35840  (-J  - 15360  (-J  ] 

(Johnson  five  — term)  (75) 

where  Ai  is  a parameter  that  relates  the  geometry  to  the  lift  coefficient. 

As  the  number  of  terms  included  in  the  series  increases,  the  loading  moves  towards  the  trailing  edge.  A 
useful  formula  for  the  finite  cavitation  number  correction  is  given  by  (Ohba  1963/1964)  as  follows. 

CL  = \J  C^)O_0  + 0.85<r2  + aCL^0  (76) 

where  Cl  is  the  required  lift  coefficient  and  Cl,<?= o is  the  lift  coefficient  for  zero  cavitation  number. 

In  Figure  20,  the  lift  to  drag  ratio,  L[D,  is  shown  for  the  NACA  4dgt  sections  designed  by  the  present 
method,  Tulin’s  two-term  sections,  and  Johnson’s  five-term  sections,  for  various  cavitation  numbers.  In  the 
same  figure  we  also  show  the  L/D  for  the  optimum  Johnson  section  as  designed  by  the  present  method.  The 
corresponding  lift  coefficients  are  shown  in  Figure  21.  For  the  design  of  Tulin  sections  and  Johnson  sections, 
the  cavitation  number  and  the  lift  coefficient  are  both  specified,  whereas  in  the  present  method,  the  cavitation 
number  and  the  lift  are  given  and  the  lift  coefficient  is  determined  as  part  of  the  solution.  A representative 
section  and  its  corresponding  pressure  distribution  is  shown  in  Figure  22. 


3.8  Application  to  3-D  elliptic  planform  hydrofoils 

As  shown  in  (Kinnas  et  al  1995),  a super-cavitating  hydrofoil  with  elliptic  planform  and  similar  sections 
along  its  span  (that  are  only  scaled  by  the  chord),  can  be  designed  by  designing  an  equivalent  2-D  section.  If  c 
and  s0  are  the  maximum  chord  and  span  of  the  elliptic  wing,  the  equivalent  2-D  lift  is  given  from: 

L-2d  = L-}l> (77) 

TTSo 

If  D2d  is  the  corresponding  minimum  value  of  the  drag  then  D3D  will  be  given  as: 


D3d  = 


Li d , ™<>D 

1/2 pU^nsl  + 4 °2D 


The  above  equation  can  be  re-written  as: 


(78) 


Djp  _ L3D  D2p 
L3D  ~ l/2pUlnsl  + L2D 


(79) 


Equation  (79)  implies  that  (L[D)3D  will  always  be  smaller  than  ( L/D)2d , he-  the  3-D  hydrofoil  will  always 
be  less  efficient  than  the  2-D  hydrofoil. 

The  actual  angle  of  attack  a3p  will  be  related  to  that  determined  from  the  2-D  equivalent  optimization 
problem  as: 


L3d 


Ot3D  = OC2D  + 


1/2  PU'^TTS20 


(80) 
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Figure  16:  Accuracy  of  the  interpolation  : Supercavitation  (Cd  vs  a),  Range  of  parameters:  a = 0.3  ~ 4.0  deg 
(6  points),  fo/c  = 0.0  ~ 0.05  (6  points),  fp/c  = 0.7  ~ 0.9  (5  points),  Ifc  = 1.1  ~ 1.3  (4  points) 


1 2 a 3 4 


Figure  17:  Sensitivity  of  the  optimum  solution  on  initial  guesses  , L = 30,000 (N/m),a  = Q.6,f7oo  = 18.2(m/s) 
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Figure  18:  L[D  for  partially  or  super-cavitating  foils  designed  by  the  present  method;  L0  = 30,00Q1V/to. 


partially  cavitating 

supercavitating 


Figure  19:  Contour  plots  of  L[D  for  partially  or  super-cavitating  foils  designed  by  the  present  method. 
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Figure  20:  Comparison  of  L/D  between  designed  sections  and  existing  sections  {L  = 30000 (N/m),zmin/c? 
7.0  x 10~9  (supercavitating) 


Figure  21:  Comparison  of  Cl  between  designed  sections  and  existing  sections  (£•  = 30000 (N/rn),zmin/c3 
7.0  x 10-5  (supercavitating) 


